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The kinematics of the outer parts of three intermediate-luminosity elliptical 
galaxies have been studied using the Planetary Nebula Spectrograph. The 
galaxies' velocity dispersion profiles are found to decline with radius; dynam- 
ical modeling of the data indicates the presence of little if any dark matter 
in these galaxies' halos. This surprising result conflicts with findings in other 
galaxy types, and poses a challenge to current galaxy formation theories. 

Over the past twenty-five years, astronomers have gone from being surprised by the exis- 
tence of dark matter to the understanding that in fact most of the Universe is dominated by 
exotic non-luminous material. In the prevailing paradigm, the gravitational influence of "cold 
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dark matter" (CDM) is crucial to the formation of structure, seeding the collapse and aggrega- 
tion of today's luminous systems. An inherent consequence of this picture is that galaxies have 
massive, extended CDM halos. Indeed, such halos are evident around spiral galaxies, where the 
rotational speeds in their extended cold gas disks do not decrease outside the visible stars — a 
gravitational signature of dark matter (|7]). 

The evidence for dark matter in elliptical galaxies is still circumstantial. Assessments of the 
total masses of individual elliptical systems have generally been confined to the very brightest 
ones, where the gravitational potential may be measured using x-ray emission ^ or strong 
gravitational lensing (3), and to nearby dwarfs, where the kinematics of individual stars offer 
a probe of the mass distribution (0). More "ordinary" elliptical galaxies, with luminosities 
close to the characteristic L* (= 2.2x10^°Lb q in _B-band solar units for a Hubble constant 
Hq = 70 km s^^ Mpc^^) Q, are more difficult to study because in general they lack a simple 
kinematical probe at the larger radii where dark matter is expected to dominate. The velocity 
distribution of the diffuse stellar light is the natural candidate (6), but studies have been limited 
by the faintness of galaxies' outer parts to radii ^ 2i?eff (where R^s is the galaxy's "effective 
radius", enclosing half its projected light). 

A powerful alternative is offered by planetary nebulae (PNe), which are detectable even in 
distant galaxies through their characteristic strong emission lines; once found, their line-of-sight 
velocities can then be readily determined by the Doppler shift in these lines. These objects have 
been used in the past as tracers of the stellar kinematics of galaxies (7), but the procedure of 
locating them using narrow-band imaging surveys and then blindly obtaining spectra at the iden- 
tified positions has proved difficult to implement efficiently on a large scale. We have therefore 
developed a specialized instrument, the Planetary Nebula Spectrograph (PN.S), specifically to 
study the kinematics of PNe in elliptical galaxies (8). The PN.S uses counter-dispersed imag- 
ing (a type of slitless spectroscopy) over a wide field to simultaneously detect and measure 
velocities for PNe using their [O III] emission at 500.7 nm. Its optimization for this purpose 
means that the PN.S is far more efficient for extragalactic PN studies than any other existing 
instrumentation. 

Observations with the PN.S on the 4.2-m William Herschel Telescope have allowed us to 
extend stellar kinematic studies to the outer parts of three intermediate-luminosity elliptical 
galaxies, NGC 821, NGC 3379, and NGC 4494 (Table 1). In each of these systems, we have 
measured ~ 100 PN velocities with uncertainties of 20 km s^^ out to radii of 4-6-Reff (Figs. 1- 
3). The line-of-sight velocities in the outer parts of all these galaxies show a clear decline in 
dispersion with radius (Fig. 4). A decrease in random velocities with radius has been indicated 
by small samples of PNe around NGC 3379 (9), but the more extensive data set presented here 
provides a more definitive measurement of this decline, and reveals that it also occurs in other 
similar galaxies. The new data are inconsistent with simple dark halo models (Fig. 4) and thus 
different from kinematical results for brighter ellipticals [e.g., (61^)]. 

More surprisingly, the velocity dispersion data follow simple models containing no dark 
matter (Fig. 4), showing the nearly Keplerian decline with radius outside 2 i?eff that such models 
predict, and suggesting that these systems are not embedded in massive dark halos. However, 
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a declining profile can also be the signature of a distribution of orbits that is primarily radial, 
because the observable line-of-sight component of velocity in such a system will decline with 
radius even if the magnitude of the velocity remains constant due to a surrounding dark halo. To 
assess the impact of this well-known degeneracy between mass and orbital anisotropy, we have 
constructed spherical Jeans models (see Appendix) with a fixed anisotropy, spanning a range 
from radially- to tangentially-biased extremes. We have fitted these models to the observed 
velocity dispersion profiles for the stars (I7illi2ll7?t and the PNe, inferring a benchmark quantity 
T B5, the i?-band mass-to-light ratio in solar units at 5 i?cfr (Table 1). Even with this extra degree 
of freedom in the modeling, the values of 5-17 are low compared with the typical T^s of ~ 20- 
40 found for bright elliptical galaxies (|2ll7?b. 

These models are not completely definitive, as there is no physical reason why the orbital 
anisotropy should not vary with radius. But we also have not used all the information contained 
in the kinematical data, fitting only binned root-mean- square values of the stellar and PN veloc- 
ities, rather than the full two-dimensional distribution of line-of-sight velocities as a function 
of projected radius. To overcome these remaining limitations, we have used an orbit library 
method (see App.) to model the dynamics of NGC 3379. This procedure takes a broad suite 
of parameterized mass models, allows for orbital anisotropics that vary arbitrarily with radius, 
and finds the best fit to the PN velocity data as well as to the stellar brightness distribution (I75l) 
and kinematics ( 72 ), including the higher-order velocity moments that quantify the shape of the 
velocity distribution. 

This procedure's only underlying assumptions are that Newtonian gravity applies, that the 
galaxy is in equilibrium, and that it is spherically symmetric. With respect to the latter assump- 
tion, the observed ellipticity of NGC 3379 varies from 0.14 to 0.30, so that its gravitational 
potential can be approximated as spherical to better than 10% accuracy. More serious is the 
possibility that this system contains a face-on disk component similar to an SO galaxy, as has 
been suggested for faint ellipticals in general (173b and for NGC 3379 in particular (15); such 
a cold component could be responsible for the apparent low velocities. However, kinematical 
modeling of the inner regions (17) and the weak rotation of the outer parts indicate that this 
is unlikely; furthermore, the galaxy's globular cluster system also shows a steeply declining 
dispersion profile (18), and this population is unlikely to contain a disk component. The model 
also relies on the PNe following the same spatial and velocity distribution as the main stellar 
component; this is not strictly true because oxygen abundance radial gradients make the PN 
detectability vary with radius, but we calculate that this will cause a mass underestimate of only 
~ 5%. Another potential complication is that background emission-line galaxies can masquer- 
ade as PNe. We have discarded several of these from the data set on the basis of their extended 
appearance, and one that is a >3-cr outlier; any remaining contaminants would cause an appar- 
ent increase in the velocity dispersion at large radii, which is not seen in the data so the sample 
does not appear to be significantly contaminated. Within these assumptions, the modeling pro- 
cess should uncover the range of possible combinations of mass and orbit distribution that could 
reproduce the galaxy's observed properties. 



3 



Using this method to analyze the data for NGC 3379 (Fig. Al) reveals that T^s is tightly 
constrained to be only 6.4 ± 0.6. Because population synthesis models indicate that the stars 
of NGC 3379 have T^* = 4-9 dTPb . and dynamical models of the inner regions where little 
dark matter is expected indicate T^^r =6-7 {20), this result suggests a "naked" galaxy with no 
dark matter cloaking its visible body. However, the mass distribution inferred by this modeling 
process does not quite match the light distribution (Fig. A2), so there appears to be a small 
amount of dark matter. A dark matter- less galaxy is still possible if the baryonic varies 
somewhat with radius, due to a gradient in the properties of the stellar population, or due to 
the presence of some undetectable baryonic component such as warm hydrogen gas. Such 
models may be plausible but are also somewhat ad hoc. In any case, the ratio of dark matter to 
total mass within 5 Rcs (9 kpc) is < 0.32, so the galaxy within this region is not dark matter 
dominated. Extrapolating the models to the virial radius (120 kpc), we find = 32 ± 13, 
which is consistent with the value of = 27 ± 5 found by an independent analysis of the 
kinematics of a cold gas ring at 1 10 kpc (21 ). 

Out of the five intermediate-luminosity ellipticals studied to date with extended kinematics 
[the three from this study as well as NGC 2434 and NGC 4697 (EH)], four have declining 
velocity dispersions. The fourth, NGC 4697, has T^s ~ 11 from Jeans modeling, so all four 
galaxies are consistent with T^s ^ 13. The apparent low dark matter content in these systems 
could be explained by pathological orbit structures (such as rapid strong anisotropy variations 
with radius), but the orbit library analysis rules out this possibility for NGC 3379, and in general 
such behavior has not been found in similar detailed modeling of elliptical galaxies. Latent disks 
could be to blame, but as with NGC 3379, the observed rotational speeds are low, and hiding 
them adequately would require coincidentally unpropitious viewing angles. Thus it seems most 
plausible that the apparent Keplerian decline in the velocity dispersions of all these systems is 
what it seems to be, and many ordinary elliptical galaxies are highly deficient in dark matter 
relative to other galaxy types — a possibility already hinted at by previous dynamical studies 

This result clashes with conventional conceptions of galaxy formation: in particular, if el- 
lipticals are built up by mergers of smaller galaxies, it is puzzling for the resulting systems to 
show little trace of their precursors' dark matter halos. More detailed comparisons with the 
predictions of the standard CDM paradigm are not yet possible, because the baryonic processes 
during galaxy formation are complex, and high-resolution ab initio simulations are so far un- 
able to reproduce an elliptical galaxy. However, estimating the effects of baryonic collapse on a 
CDM halo (26) using an adiabatic approximation yields T^s ~ 20. Low resolution simulations 
including baryonic processes (27), extrapolated inward in radius, also predict T55 ~ 20. These 
values conflict with those derived for the observed galaxies above. On the other hand, CDM 
predictions at the virial radius are generally supported by statistical studies of L* ellipticals us- 
ing weak gravitational lensing (I2S'I) . This apparent disparity could be resolved if the dark matter 
content of ellipticals varies, or if these galaxies have large amounts of dark matter spread out 
to still larger radii than the PNe can probe. The former scenario begs a suitable mechanism for 
dark matter depletion, and the latter violates CDM predictions for high central concentrations 
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of dark halos — a problem that is also increasingly evident in other stellar systems (EPb . 

It is apparent that some important physics is still missing from the recipes for galaxy forma- 
tion. One obvious candidate is that substantial portions of these galaxies' dark matter halos have 
been shed through interactions with other galaxies. Such stripping has been inferred for ellipti- 
cals near the centers of dense galaxy clusters but the galaxies we studied here are in much 
sparser environments, where substantial stripping is not expected to have been an important 
process. Crucial to understanding the incidence and origin of this low dark matter phenomenon 
will be results for a large sample of ellipticals with a broad range of properties, including dif- 
fering environmental densities, which could be a key factor in determining halo outcomes; the 
continuing PN.S observing program will provide this sample. 
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Galaxy 


Type 


D (Mpc) 


-Reff 
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Observations 


A^(PNe) 
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E2 


24 


50" 


-20.5 


Sep 01, llh 


104 


13-17 


NGC 3379 


El 


11 


35" 


-20.0 


Mar 02, 3h 


109 


5-8 


NGC 4494 


EO 


17 


49" 


-20.6 


Mar 02, 4h 


73 


5-7 



Table 1: Properties of three galaxies observed with the Planetary Nebula Spectrograph. 
Columns show galaxy name, type, distance (ETl), effective radius (32), absolute magnitude, 
date and effective integration time of observations (normalized to I" seeing), number of PNe 
detected, and the range of 5-band mass-to-light ratio at 5 effective radii (in solar units) deter- 
mined from anisotropic Jeans models. The mass-to-light ratios are systematically sensitive to 
the assumed distance — e.g., if D = 33 Mpc for NGC 821 {33), T^s is decreased by 30%; and 
if D = 11-13 Mpc for NGC 4494 Q3.M}, is increased by 30-60%. The characteristic 
magnitude is = -20.4 ©. 
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Fig. 1. NGC 3379 with 109 planetary nebula line-of-sight velocities relative to the systemic 
velocity, as measured with the William Herschel 4.2-m telescope and the PN.S instrument. Red 
crosses represent receding velocities, and blue boxes are approaching, where the symbol sizes 
are proportional to the velocity magnitudes. The background image is from the Digitized Sky 
Survey (35), and the field of view is 8.'4x8.'4 = 26x26 kpc = 14x 14 Rcs- 
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Fig. 2. NGC 821 with 104 planetary nebula velocities. See Fig. 1 for further description. The 
field of view is lO'x 10'= 70x70 kpc = 12xl2i?efr- 



9 




1 arcmln; 



Fig. 3. NGC 4494 with 73 planetary nebula velocities. See Fig. 1 for further description. The 
field of view is 10'xlO'= 50x50 kpc= 12xl2i?efr- 
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Fig. 4. Line-of-sight velocity dispersion profiles for three elliptical galaxies, as a function of 
projected radius in units of the effective radius. Open points show planetary nebula data (from 
PN.S); solid points show diffuse stellar data (I77]|72]|221l- The vertical error bars show 1 a uncer- 
tainties in the dispersion, and the horizontal error bars show the radial range covered by 68% of 
the points in each bin. Predictions of simple isotropic models are also shown for comparison: a 
singular isothermal halo (dashed lines) and a constant mass-to-light ratio galaxy (dotted lines). 
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Appendix: Methods 



Jeans models 

The non-rotating spherical Jeans equations link the second moments o"j (r) of a system's velocity 
distribution to its mass M (r) and luminosity z/(r) distributions (OHb . We use the formulation 



where and are the radial and tangential velocity dispersions, and (3 = l — ag/a^ quantifies 
the (unknown) velocity anisotropy, We take the luminosity distribution to be the Hemquist ( 1771) 
profile z/(r) = L{2nr)^^{r + a)~'^ corresponding to the known effective radius of the galaxy 
(i?eff = 1.8153a), and parametrize the dispersion profile by the form 



(thus allowing for a dark halo with a range of radial mass distributions), which is projected and 
fitted to the dispersion data while optimizing for the parameters {vq, tq, a, S}. The anisotropy 
is set to be constant and is varied between the reasonable extremes of /? = —0.5 (tangential) to 
P = 0.5 (radial). We then solve the Jeans equation for the mass profile M(r). 

These models are suited to exploring the range of outer galaxy mass profiles, and not for 
reproducing the fine structure of the dispersion data inside 1 R^s. They are thus illustrative 
of the range of T b5 associated with the anisotropy uncertainty, and these are the T b5 values 
reported in Table 1 ; the statistical uncertainty is of the same order of magnitude as this system- 
atic uncertainty. As discussed for NGC 3379, the assumption of spherical symmetry should be 
accurate enough for mass estimates at the 10% level. 

Orbit models 

Our orbit library method is documented in (38). In summary, it begins with an assumed form 
for the mass distribution, and in its gravitational potential calculates a library of orbits, of all 
types from radial to circular. These individual orbits are then combined non-parametrically to 
find the best fit to the observational constraints. Regularization is used to stabilize the solution, 
but at such a low level as to avoid biasing them. The fit to the data is quantified by a likelihood 
function £, and 68% confidence limits in model parameters are given by A In £ = 1/2. Because 
the degrees of freedom are ill-defined with this method, the absolute value of In C is not very 
informative, but the diffuse stellar component of the fits (photometric and kinematic) has a 
typical = 8-9 for 129 data points, and the models qualitatively reproduce the data properties. 




(1) 




(2) 
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We use a grid of models consisting of a constant mass-to-light ratio Hemquist model galaxy 
(ITTb plus a CDM simulated dark halo (l23t . with the mass density distribution: 



27rGr(r + a)3 ATcGr{rs + ry 

The mass model parameters are thus the halo core radius r^, the mass of the stars v^, and the 
relative masses of the stars and halo vl/vl. Our grid covers values of from 113" to 700" (6 
to 36 kpc), from 575 to 625 km s^^, and vl/vl from to 1.1. 

The best-fit solution for NGC 3379 is shown in Fig. Al, and the range of permitted solutions 
is illustrated by Fig. A2. We find T^g (in solar units) to be 6.4±0.6 at 185" (9 kpc), 7.9± 1.1 at 
270" (14 kpc), and 32 ± 13 at 2460" (120 kpc). Allowable solutions within the 68% confidence 
limits include {r„ v,, vl/vl} = {199", 603, 0.15}, {280", 610, 0.15}, {350", 610, 0.12}, {350", 
603, 0.20}, {500", 610, 0.12}, {500", 603, 0.20}, and {500", 603, 0.30}, where v, is in km s'^. 
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Fig. Al. Line-of-sight velocity distribution with radius in NGC 3379. The error bars show stel- 
lar dispersion data (72); the points show the PN velocity data. The background shading shows 
the projection of the best-fit model, multiplied by the stellar luminosity at each radius; the dotted 
line shows its dispersion. 
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Fig. A2. Intrinsic circular velocity profile with radius, fc(r) = y GM{r)/r, of NGC 3379. The 
solid lines and shaded area show the region permitted by orbit modeling. The dotted lines show 
models which are ruled out at the l-a level: the bottom one shows a constant mass-to-light ratio 
(M/L) model, the upper three show models with more dominant dark halos. Each model is 
labeled with its fi-band mass-to-light ratio (in solar units) at 5 effective radii. Thus we see that 
the shape of fc(r) for a constant M/L galaxy falls slightly outside the allowed range. 
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